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ABSTRACT 

We study ~330 massive (M* > 10 95 M Q ), newborn spheroidal galaxies (SGs) around 
the epoch of peak star formation (1 < z < 3), to explore the high-redshift origin of SGs 
and gain insight into when and how the old stellar populations that dominate today's 
Universe formed. The sample is drawn from the HST/WFC3 Early- Release Science 
programme, which provides deep 10-filter (0.2 - 1.7 /im) HST imaging over a third 
of the GOODS-South field. We find that the star formation episodes that built our 
SGs likely peaked in the redshift range 2 < z < 5 (with a median of z ~ 3) and have 
decay timescales shorter than ~1.5 Gyr. Starburst timescales and ages show no trend 
with stellar mass in the range 10 9 ' 5 M < M* < 10 10 ' 5 M Q . However, the timescales 
show increased scatter towards lower values (< 0.3 Gyr) for M* > 10 10 5 M©, and 
an age trend becomes evident in this mass regime: SGs with M* > 10 11 " 5 M are 
~2 Gyrs older than their counterparts with M* < 10 10 5 M Q . Nevertheless, a smooth 
downsizing trend with galaxy mass is not observed, and the large scatter in starburst 
ages indicates that SGs are not a particularly coeval population. Around half of the 
blue SGs appear not to drive their star formation via major mergers, and those that 
have experienced a recent major merger, show only modest enhancements (~40%) in 
their specific star formation rates. Our empirical study indicates that processes other 
than major mergers (e.g. violent disk instability driven by cold streams and/or minor 
mergers) likely play a dominant role in building SGs, and creating a significant fraction 
of the old stellar populations that dominate today's Universe. 

Key words: galaxiesdormation - galaxies: evolution - galaxies: high-redshift - galax- 
ies: elliptical and lenticular, cD - galaxies: interactions 



1 INTRODUCTION 

Massive galaxies are central to our understanding of the vis- 
ible Universe. Locally, the massive galaxy census is dom- 
inated by systems with spheroidal morphology. Hosting 
more than ^50% of the stellar mass at the present day 
(e.g. iBernardi et"alll2003h , these massive spheroidal galax- 



ies (SGs) are key laboratories for studying the evolution 
of g alaxies over c osmic time. The i r red optical colours 
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(e.g. IJorgensen et al.lll99fj; iFranxl Fl995l ; ISaglia et al J 1 19971 ; 
iForbes et all 1 19981 ; iPeeblej I2002T ) and chemical evidence 
for short star-forma tion timescales (e.g. iTrager et all l2000l : 
iThomas et "all 120051 ). indicate that the bulk of their stellar 
mass formed rapidly at high redshift (z > 1). Neverthe- 
less, the young (luminosity-weighted) ages observed in some 
SGs are a signature of recent (< 5 Gyr) s tar formation (e.g . 
Ellis et al.ll200ll; iMenanteau et all 1200 ll; iNelan et al.ll2005l ; 



Schiavorj2007l ; lTrager et alj200g| : lGraves et al.ll2009T ). which 
can be accurately quantified using a sensitive probe of young 
stars, such as the rest-frame ultraviolet (UV; 1200-3000 A). 
Recent UV studies have shown that, while old stars do dom- 
inate today's SGs, a significant minority (~20%) of the stel- 
lar mass in these galaxies is formed at z < 1 (Kaviraj et al. 
2007a; 2008, see also Ferreras & Silk 2000a; Yi et al. 2005; 
Salim & Rich 2010; de la Rosa et al. 2011; Rutkowski et al. 
2012), via minor mergers between SGs and gas-rich dwarfs 
(Kaviraj et al. 2009a; Kaviraj et al. 2011, see also Tal et al. 
2012; Newman et al. 2011). 

While it is clear that the dominant stellar populations 
in SGs are old, remarkably little is known about exactly 
when and how these stars formed in the early Universe. 
The key questions are: (1) At what redshifts and over what 
timescales did they form? (2) What were the principal mech- 
anisms that drove this star formation? While the former 
question has been explored using spectro-photometric stud- 
ies of local SGs, it is challenging to construct a detailed star 
formation history for the old stars using local SGs, because 
the available spectro-photometric indicators are insensitive 
to old stellar populations of different ages. Optic al colo urs, 
for example, evolve slowly after 4-5 Gyr (e.g. lYil 12003! ) . A 
further complication is that all indicators are affected by 
recent star formation, typically by an uncertain amount in 
any given galaxy. While old stellar populations do domi- 
nate SGs by mass, measured spectro-photometric quantities 
are luminosity-weighted. Since the luminosity-weighting of 
young stars is higher than their old counterparts, even small 
amounts of recent star formation can have a disproportion- 
ately large impact on the indicator in question. Thus, while 
it is possible to put lower limits on the age of the old stellar 
populations in today's SGs, the exact shape of the star for- 
mation history of these old stars remains elusive. It is worth 
noting that the average timescale of star formation in indi- 
vidual galaxies can be estimated with reasonable precision 
using [a/Fe] ratiofl Indeed, the super-solar [a/Fe] ratios 
observed in mos t SGs suggest star-formation timescales less 
than ~ 1 Gyr dTrager et all |2000| ; iFerreras fc Silk! l2000bl ; 
IThomas et al.ll2005l ) 

The latter question (2) has been the focus of 
some debate in the literature. While the classical view 
has been that SGs are products of 'major' (roughly 
equal-mass) mergers betwe en mass i ve sp ir al galaxies in 



the early Universe (e.g. iToomre 19771; White! 
iBarnes fc Hernquistl 1 19921 ; ISomerville fc Primackl 



1978 



1!)!)!) 



1 a elements are primarily provided by prompt Type II super- 
novae, which explode almost instantaneously (~ 10 6 yr) on cos- 
mological timescales. Fe, on the other hand, is provided by Type 
la supernovae (S N), which emerge aft er typical time delays of less 
than a Gyr (e.g. IThomas et al.lll999l) ■ The [a/Fe] ratio therefore 
indicates the ratio of the overall timescale of star formation in a 
galaxy to the onset timescale of the Type la SN. 
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Figure 1. Redshift distribution of the SGs in this study. 



Relaxed spheroids 



162 (z= 1.54) 




Figure 2. Examples of relaxed SGs in our dataset. The redshift 
of each galaxy is indicated above the image, next to its ID. We 
show both the YJH colour image and its negative. 



ICole et all l2000l : lHatton et ail 120031 : ISpringel et all 120051) 

recent hydr o dyna m ical simulations ( e.g. | Dekel et al. 120091 ; 
iKeres et all |2009| ; iDevriendt et alj |2010| ). coupled with 
the observed paucity of major merger s at high redshift 
(e.g. iGenzel et al.ll2008l ; lLaw et al.ll2009l ; Ijogee et al.ll2009l ; 
iKavirai et al. 2011 ). suggest that the progenitors of SGs 
may not predominantly be major mergers but rather clumpy 
disks, in which star f ormation is fed by cold streams and 
minor mergers (see e. g. lBirnboim fc Dekellliool; IKeres et all 
120051 ; iDekel et al.ll2009l ; ICeverino et al.ll2010l . l2012| ), with the 
clumps eventually coalescing to form a spheroid. 

Given these significant open issues, a detailed empir- 
ical understanding of the formation of the dominant, old 
stars in today's SGs is highly desirable, both as a route 
to understanding the drivers of star formation in the early 
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Universe, and as an important test of our emerging theo- 
retical models for the high-redshift Universe. Given the lim- 
ited utility of local SGs in studying these issues, an ideal 
method for probing questions (1) and (2) above is to study 
newborn SGs around t he epoch of peak cosmic star for- 
mation (1 < z < 3, see Madau et al.lll998l ; I Hopkins! 12001 
IHopkins fc Beacoml lioog T Quantifying the star formation 
histories of these systems then allows us to directly probe 
the formation of the old stars in today's SGs, free from the 
contamination of the galaxy spectrum by intervening star 
formation episodes at late epochs. The underlying assump- 
tion is that newborn SGs at high redshift are the ances- 
tors of their present-day counterparts, i.e. that SGs do not 
routinely transform into disks over cosmic time. Both the- 
oretical and observational work on the evolving morpholog- 
ical mix of the Universe suggests that this is indeed a rea- 
sonable assumption to make (e.g. Butcher fc Oemlei 19841 ; 
iBaugh et al1ll99rj ; ICouch et al.lll998l ; ISmith et al.ll2005l ). 

While the analysis of newborn SGs at z > 1 is a com- 
pelling pr oject, past studies o f statistically large samples of 
SGs (e.g. lFerreras et al.ll2009l ) have typically fallen short of 
this epoch, because deep, survey-scale HST imaging in the 
near-infrared (NIR), required for obtaining rest-frame opti- 
cal galaxy morphologies at z > 1, has, until recently, been 
lacking. However, fol lowing the GOODS-NICMOS survey 
(jConselice et al.ll201lh . a new generation of NIR surveys us- 
ing the HST's WFC 3 - e.g. the WFC3 Earl y-Release Science 
(ERS) programme (|Windhorst et al.ll201lh an d the ongoing 
CAN DELS campaign ( Koekemoer et alJ 20 111 : iGrogin et al.l 
l201ll ) - are providing unprecedented large-scale access to 
rest-frame optical galaxy data at z > 1, making them ideal 
datasets for the study of newborn SGs at high redshift. Here, 
we employ WFC3 ERS data and study SGs around the 
epoch of peak cosmic star formation, to explore the high- 
redshift origin of these galaxies and gain insight into when 
and how the old stellar populations that dominate today's 
Universe formed. We leverage the unprecedentedly deep NIR 
imaging to morphologically identify SGs and explore the role 
of major mergers in their evolution, and exploit the uniquely 
wide ERS wavelength baseline (0.2 - 1.7 (im) to quantify 
galaxy star formation histories via the available rest-frame 
UV-optical photometry. 

This paper is organised as follows. In Section 2, we 
briefly describe the ERS dataset and the morphological se- 
lection of SGs using visual inspection. In Section 3, we de- 
scribe the calculation of star formation histories and rest- 
frame photometry for individual SGs in our sample. In Sec- 
tion 4, we study the photometric properties of the new- 
born SGs, use their derived star formation histories to con- 
strain the formation epoch of their stellar mass and explore 
the mechanisms by which this stellar mass is likely being 
formed. We summarise our findings and outline avenues for 
future work in Section 5. Throughou t this paper we use th e 
WMAP7 cosmological parameters (|Komatsu et al.l l201ll). 
All p hotometry is presented in the AB system ( Oke fc Gunnl 
Il983ft . 



GALAXY SAMPLE AND 
MORPHOLOGICAL CLASSIFICATION 



Disturbed spheroids 



1296 (z= 1.65) 




Figure 3. Examples of disturbed SGs in our dataset. The redshift 
of each galaxy is indicated above the image, next to its ID. We 
show both the YJH colour image and its negative. 



The WFC3 ERS programme has imaged around one-third 
of the GOODS-South field with both the UVIS and IR chan- 
nels of the HST/WFC3. The observations, data reduction, 
and instrument performance are described in detail in Wind- 
horst et al. (2011), and highlighted here. The goal of this 
part of the ERS programme was to demonstrate the sci- 
ence capabilities of the WFC3 for studying intermediate- 
and high-redshift galaxies in both the UV and NIR, by ob- 
serving a portion of t he well-studied GOODS-South field 
ijGiavalisco et al.l 120041 ). The field was observed for 104 or- 
bits, with 40 orbits in the UVIS channel and 60 orbits in 
the NIR channel. The UVIS data covered ~55 arcmin 2 , in 
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K z < 1.5 1.5 < z < 3 

fsG N DS g/N sg fsG N DSG /N SG 

Log M« > 9.5 0.35 ±0 - 02 0.42± - 04 0.32± 002 0.35± 04 

9.5 <Log M, < 11 0.36± 03 0.43 ±0 05 0.36± 002 0.36± 04 
Log M* > 11 0.56± 011 0.33± 012 0.58± 013 0.18± 08 

Table 1. The fraction of galaxies that are SGs (fso) m various stellar mass (units of Mg) and redshift ranges, and the fraction within 
the SGs that are disturbed (Ndsg/Nsg)- For the most massive galaxies (lO n M0-lO 12 M ), for example, the fraction of SGs is ~60%, 
similar to the value found by Cameron et al. 2011, who also employed visual inspection for identifying ETGs. The total number of 
galaxies in this study is 818 (736 in the range 9.5 <Log M* < 11 and 82 with Log M, > 11). 



each of the F225W, F275W and F336W filters, with relative 
exposure times of 2:2:1. The IR data covered ~45 arcmin 2 
using the F098M (Y s ), F125W (J), and F160W (H) filters 
with equal exposure times of 2 orbits per filter. The data 
were astrometrically aligned with a version of the GOODS- 
S HST/ACS data (v2.C0; Giavalisco et al. 2004) that was 
rebinned to have a pixel scale of 0'.'09 per pixel. Together, 
the data provide 10-band HST panchromatic coverage over 
0.2 - 1.7 fim, with a 5cr point-source depth of AB ~ 26.4 
mag in the UV and AB ~ 27.5 mag in the IR. 

Our focus in this paper is the subset of 818 
ERS galaxies that have either spectroscopic or photo- 
metric redshifts in the range 1 < z < 3. Photo- 
metric redshifts are calculated using the EAZY code 
ijBrammer et alJ I2008T I on the 10-band WFC3/ACS pho- 
tometric catalogue (Cohen et al. in prep). Spectroscopic 
redshifts are drawn from the liter ature, from spectra 
taken using the Ver y Large Telescope (iLe Fevre et al.ll2004l: 



Szokolv et all 12004 iMignoli et all 120051: iRavikumar et all 
20071 : IVanzella et all [20081: iPopesso et all 120091 ) . the Keck 



telescopes (IStrolger et all 120041) and the HST ACS grism 
jDaddi et all 120051 : iPasauali et al l I2006J : iFerreras et all 
2009). For the analysis that follows (e.g. in Section 3), spec- 
troscopic redshifts are always used where available. Figure 
[1] shows the redshift distribution of the SG sample in this 
study. 

Each galaxy in this sample is morphologically classified 
via visual inspection of its WFC3/NIR images, which trace 
the rest-frame optical wavelengths at z > 1. While morpho- 
logical proxies, such as concentration, asymmetry, dumpi- 
ness, M20 and the Gini coefficient have b een widely used 
to classify galaxies in large surveys (see e.g. Abraham et alJ 

20041 ). the per^ 



1 19961 . 120031 : IConselice et aljlifjoll : lLotz et al 



formance of these met hods is typically cali brated against 
visual inspection (e.g. I Abraham et al.l 1 19961 ). which offers 
better pr ecision and reliability in morphological classifica- 
tion (e.g. lLiskerll2008l : iRobaina et~afll2009l ; iKartaltepe et al.1 
l201Ct ). Past HST studies of the high-redshift Universe 
(including recent efforts using the WFC3) have com- 



phologies (e.g. 


Elmegreen et al. 20051: Ferreras et al.l 


2005; 


Bundv et all 2005: ICassata et al. 120051: Jogee et al.l 


2009; 


Robaina et al.1 


2009: iKavirai et al.l 201 ll ICameron et al. 


201ll; Kocevski et alj|2012f), using rest-frame optical imag- 



ing that has similar or fainter surface-brightness limits com- 
pared to the images employed here. 

In this paper, we visually inspect YJH composite im- 



http: / / archive, stsci . edu /pub /hlsp /goods/v2/ 



Reloxed SGs 



> 4 ■ 
I 

> 



• — s V .n .Vf ••• 



-0.5 



0.0 0.5 
(B-V). 



1.0 



> 4 ■ 
I 

> 



Disturbed SGs 



I 



l 12 
11.3 

10.7 

10.1 

9.50 



-0.5 



0.0 



0.5 
(B-V). 



1.0 



Figure 4. Rest-frame UV/optical colours of the ERS galaxies in 
the redshift range 1 < z < 3. The top panel shows SGs that are 
relaxed i.e. do not show any morphological disturbances, while 
the bottom panel shows SGs that are morphologically disturbed. 
The filled coloured circles represent SGs, while the small black 
dots represent late-type galaxies. The colour coding indicates the 
stellar mass of the SG (see legend). The errors in the stellar masses 
are typically better than 0.3 dex. 



ages, scaled using the asinh method o f lLupt on et al.l l|2004h . 
to determine (1) the morphology of the galaxy and (2) 
whether it exhibits morphological disturbances, i.e. tidal 
features indicative of a recent merger. Using YJH com- 
posites (instead of monochrome images in one filter) max- 
imises the rest-frame optical information in the image, and 
facilitates the identification of tidal features. While previ- 
ous studies have successfully performed visual inspection 
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Figure 5. The offset between the rest-frame colours of individual 
SGs and that of a dustless, solar-metallicity instantaneous burst 
at z = 5. The dotted line shows how the offset will change if 
we consider a more realistic instantaneous burst that is folded 
with the median value for the internal extinction of the SGs 
(Eb—v ~ 0.2 mag). The arrow indicates the reddening vector for 
this extinction value. Black and red circles indicate relaxed and 
disturbed SGs respectively and small black dots indicate late- type 
galaxies. This comparison, intended only as a guide, indicates that 
none of the SGs in our sample are likely to have completed their 
stellar mass assembly by z ~ 5. 



of sim ilar datasets to H ~ 25 mag (e.g. I Cameron et al.l 
l201ll >. we restrict our galaxy sample to a more conservative 
H ~ 24.5 mag (in this magnitude range the galaxy sample 
with M* > 10 9 ' 5 Mq is expected to be complete; Windhorst 
et al. 2011). Galaxies are classified into two morphological 
classes: SGs and a (broad) class of 'late-types', which in- 
cludes all objects that are not spheroidal, e.g. disks, clumpy 
galaxies, irregulars, etc. Here, we focus solely on the SGs, 
splitting them into 'relaxed' systems, which do not exhibit 
morphological disturbances indicative of recent interactions, 
and 'disturbed' systems, which do (see Figures [2] and [3] for 
examples of relaxed and disturbed SGs respectively). The 
visual classification process identifies 328 SGs in our sam- 
ple. As we describe in Section 4 below, the disturbed SGs 
are likely to have experienced recent major mergers. 

For the analysis that follows, we study ERS galaxies 
that have redshifts in the range 1 < z < 3, ff-band magni- 
tudes brighter than 24.5 (where morphologies are reliable) 
and stellar masses greater than 10 9 ' 5 Mq (the derivation of 
stellar masses is described in Section 3 below) . Table 1 sum- 
marises the fraction of galaxies that are SGs (fsG) in various 
stellar mass and redshift ranges and the fraction within the 
SGs that are disturbed (Ndsg/Nsg)- For example, for the 
most massive galaxies (M* > 10 11 Mq), the fraction of SGs 
is ~60%, similar to the value found by Cameron et al. (2011), 
who also used visual inspection of WFC3 images for mor- 
phological classification. In the mass range considered here 
(M* > 10 9 - 5 M Q ), around a third of SGs exhibit morpholog- 
ical disturbances. The disturbed fraction decreases with in- 
creasing redshift, plausibly due to greater cosmological dim- 
ming of the tidal features at larger distances. 



3 STAR FORMATION HISTORIES AND 
REST-FRAME PHOTOMETRY 

In this section, we describe the estimation of star forma- 
tion histories (SFHs) and rest-frame photometry that are 
used in Section 4. We compare the observed photometry 
of each galaxy to a large library of synthetic photometry. 
The synthetic library is constructed using exponentially- 
decaying model SFHs, described by five free parameters: the 
age (T), e- folding timescale (r), metallicity (Z) and inter- 
nal extinction (Eb-v) of the starburst and the total stellar 
mass formed (M*). We vary T between 0.05 Gyrs and the 
look-back time to z = 20 in the rest-frame of the galaxy, r 
between 0.01 Gyrs (approximately an instantaneous burst) 
and 9 Gyrs (approximately constant star formation), Z be- 
tween 0.05 Zq and 2.5 Zq and Eb-v between and 2 
mag. Synthetic photometry is generated for each model SFH 
by folding it with the stellar mo dels of Chariot fc Bruzua l 
(2007) (an updated version of the lBruzual fc Chariot! (|2003h 
models), through the correct WFC3 and AC S filter through- 
puts. The empirical law of ICalzetti et all l|2000h is used 
to compute the dust-extincted spectral energy distributions 
(SEDs). The synthetic library has ~ 2.5 million individual 
models. Since our galaxies span a large range in redshift 
(1 < z < 3), equivalent libraries are constructed at redshift 
intervals of Sz = 0.01. 

For each galaxy, the values of the free parameters (T, 
t, M„, Z, Eb-v) are estimated by comparing its observed 
photometry to every model in the synthetic library that 
is closest to it in redshift. The likelihood of each model, 
exp(— X 2 /2), is calculated using the value of x 2 i computed 
in the standard way. From the joint probability distribu- 
tion (which is a function of all five free parameters), each 
individual parameter is marginalised^ to extract its one- 
dimensional probability density function (PDF). The me- 
dian of this PDF is taken as the best estimate of the param- 
eter in question and the 25 and 75 percentile values (which 
enclose 50% of the probability) are used to calculate an un- 
certainty on this estimate. The K-corrections required to 
construct rest-frame photometry for each galaxy are calcu- 
lated using the best-fit model SED (i.e. where the value of 
X 2 is a minimum). 

It is worth noting that the accuracy of the derived star 
formation histories is aided by the availability of photometry 
covering both the rest-frame UV and optical wavelengths. 
The rest-frame UV, in particular, offers almost an order of 
magnitude greater sensitivity to recent star formation than 
the optical wavelengths (Kaviraj et al. 2007a), and is largely 
free from the effects of the age-metallicity degeneracy (Kavi- 
raj et al. 2007b). The rest-frame UV colours are therefore a 
strong indicator of how truly quiescent a galaxy is. While the 
rest-frame UV photometry constrains the gradual quench- 
ing/decline of star formation in our SGs, the rest-frame op- 
tical constrains the epoch of stellar mass assembly. 



3 The one-dimensional marginalised probability distribution 
(PD) for age (T), for instance, is obtained from the joint PD by 
integrating out the other parameters. If P(T, r, M, Z, Eb_v|D) is 
the joint PD (given the data D), then the marginalised PD for T, 
P(T|D) = //// oo P(T,r,M»,Z,E B .v|D) dr dM* dZ dE B _v 
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Figure 6. Typical remnants of mergers, drawn from a hydrodynamical cosmological simulation, with mass ratios of 1:2 (left), 1:5 
(middle) and 1:10 (right). The colour-coding indicates the surface-brightness (mag arcsec - 2 ) in the H-band (the surface brightness 
limit of our WFC3/IR images is ~ 26 mag arcsec -2 ). To allow comparison with the blue SGs in our sample, the remnants are 
'observed' while they are in the UV-optical blue cloud, typically ~0.3-0.4 Gyr after the merger is complete. The simulations suggest 
that tidal debris around remnants of mergers that have mass ratios less than ~1:5 are very unlikely to be observable at the depth of 
the ERS images across the redshift range considered in this study. In other words, blue SGs that do exhibit such tidal features arc 
likely to be major-merger remnants (see text for more details). 



4 THE MASS ASSEMBLY OF NEWBORN 
SPHEROIDS 

We begin by exploring the photometric properties of our 
ERS dataset. Figure [4] presents the rest-frame UV-optical 
colours of the ERS galaxies. T he rest-frame NUV wave- 
lengths, based on the GALEX l|Martin et al.l 120051 ) NUV 
filter, are centred at ~2300 A. Relaxed and disturbed 
SGs span the entire colour space occupied by the galaxy 
population, with ~40% and ~60% of relaxed and dis- 
turbed SGs respectively populating the 'blue cloud' (B — 
V < 0.5 and NUV - V < 2.5). This is consistent with 
the large range in rest-frame (UVJ) optical colours ob- 
serve d in the massive gala xy population at these redshifts 
(e.g. IWhitaker et al.l |2010| ) and the gradual, rather than 
abrupt, decline of star formation implied in most massive 
objec ts in this redshift r ange by their spectral features 
(e.g. [ Cimatti et al.l 120081: Ivan Dokkum fc Brammej |2010| ; 



iKriek et al.ll201ll : Ivan Dokkum et al-I^OHl ). In Figure [5] we 
show the offset between the rest-frame colour of individual 
SGs and that of a dustless, solar-metallicity instantaneous 
burst at z — 5. This comparison, intended only as a guide, 
indicates that very few of the SGs are consistent with having 
completed their star formation by z ~ 5. Note that, adding 
the median internal extinction in SGs derived in Section 3 
(Eb-v ~ 0.2 mag, see reddening vector in Figure [5]) to this 
instantaneous burst, only reinforces this conclusion. The star 
formation in SGs at 1 < z < 3 is, therefore, either ongoing 
or recently completed. 

It is worth noting that around half of the SGs that 
inhabit the blue cloud appear relaxed. Given that tidal fea- 
tures are readily visible in a significant fraction of objects, 
we explore whether the lack of such features in the blue re- 
laxed SGs implies that the star formation in these systems is 
not being driven by a recent major merger. Indeed, tidal de- 
bris from a recent merger is most readily visible in the early 
stages of relaxation, when star form ation remains strong an d 
the galaxy is in the blue cloud (e.g. ICarpineti et al.ll2012T ). 

To explore the detectability of tidal features in our 
high-redshift SGs, we appeal to merger remnants of vari- 



1 .0 



0.6 



0.4 



0.2 



0.0 



Reloxed SG (bluest 30%) 
Disturbed SG (bluest 30%) 



-9.4 -9.2 -9.0 -8.8 -8.6 -8.4 
Log sSFR [yr"'] 



-8.2 



3.0 



Figure 7. Specific star formation rates (sSFRs) in the bluest 
30% of relaxed (black) and disturbed (red) SGs (in the NUV - r 
colour). Median values arc shown using the dotted lines. The gen- 
eral enhancement in sSFR due to a merger appears to be modest, 
with the median sSFR of the bluest 30% of disturbed SGs being 
around 0.15 dex (~40%) higher than the corresponding median 
in their relaxed counterparts. Median values are shown by the 
dotted vertical lines. Note that the SFR is an average value over 
the last 10 7 yr. 



ous mass ratios drawn from a hydrodynamical cosmolog- 
ical si mulation . The simulation, which is described in de- 
tail i nlPeiranil (I2010f) . was performed using the GADGET2 
code (|Springe]H2005T ) with added prescriptions for star for- 
mation, feedback from Type la and II supernovae, a UV 
background and metal enrichment. The dark matter and 
baryonic particle resolutions are tudm = 7.4 x 10 6 /i _1 Mq 
and m gas = m star = 1.5 x 10 6 fo -1 M fl respectively. We refer 
readers to Section 2.1 of lPeiranH l|2010h for further details of 
the simulation. 

To compare with the blue SGs in our sample, we explore 
the surface brightness of tidal features in merger remnants 
at z ~ 1.25, the midpoint of the lower redshift bin in Table 
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Figure 8. TOP LEFT: Derived starburst age (T) vs. stellar mass in individual SGs. The yellow shaded region in the top-left panel 
shows a progressive 'one-sigma' fit to the data. The red line indicates the mean value, while the yellow shaded area indicates the 
region that encloses 68% of the galaxies. The typical errors in the stellar masses are better than 0.3 dex. TOP RIGHT: Derived 
starburst timescale (r) vs stellar mass in individual SGs. BOTTOM LEFT: Distribution of stellar mctallicities in individual SGs. 
BOTTOM RIGHT: Distribution of internal extinction in individual SGs. Median values arc shown using the dotted vertical lines. 
The error bars in the panels for metallicity and internal extinction indicate average uncertainties in the derived parameters across 
our SGs. 



1. We focus on the lower redshift bin in order to minimise 
the impact of cosmological dimming - recall that the dis- 
turbed SG fraction decreases with redshift in Table 1, plau- 
sibly due to cosmological dimming of the tidal features with 
increasing redshift. In Figure [6] we present typical examples 
of remnants of mergers that have ratios of 1:2, 1:5 and 1:10. 
The merger remnants are 'observed' while they are in the 
UV-optical blue cloud (to ensure consistency with our blue 
observed SGs) and the synthetic images are constructed by 
combining the simulation o utputs with the cor rect WFC3 
filter throughputs, following iPeirani et al.l |2010T ). Note that 
our aim here is to study the surface brightness of tidal fea- 
tures around merger remnants and not a detailed explo- 
ration of the morphology of the remnants themselves, which 
would benefit from a simulation with higher resolution. 

We find that, while major mergers produce strong fea- 
tures which will be readily visible, mergers with mass ra- 
tios of 1:5 are only marginally detectable at z ~ 1.25, 
given the surface-brightness limit of extended objects in our 
WFC3/IR images (~ 26 mag arcsec -2 ). Mergers with mass 



ratios lower than this value (see e.g 1:10 example in Fig- 
ure [6} will produce tidal features that will not be visible in 
our images. At the upper limit of our redshift range (z~3), 
only mergers with mass ratio ~1:2 will be visible. In sum- 
mary, it is likely that only tidal features produced by major- 
mergers are likely to be visible across the redshift range con- 
sidered in this study. In other words, if a blue SG in our high- 
redshift sample has experienced a recent major merger, then 
the tidal debris should be visible in the ERS images. The 
dominant star formation mechanism in the blue relaxed SGs 
therefore appears unlikely to be a major merger, implying 
that these systems are driving their star formation via other 
mechanisms, such as minor mergers, or are the remnants 
of recently collapsed (cold-stream-fed) cl umpy disks, as en- 
visaged in current the oretical work (e.g. iDekel et alj|2009l : 
lElmegreen et al]|2008f ). Note that minor-merger-driven star 
formation could, in principle, be distinguished from disk col- 
lapse by the presence of tidal features. However, as Figure [6] 
indicates, our WFC3 images are too shallow to detect tidal 
debris from minor mergers. While it is challenging to disen- 
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tangle these two processes using current data, it might be 
possible to perform this exercise using deeper data, either 
using the WFC3 itself or from future instruments such as 
the JWST or the extremely large telescopes. 

It is worth noting further that the relaxed SGs in the 
blue cloud have a similar UV-optical colour distribution in 
Figure [4] as their disturbed counterparts. Since colours re- 
flect the (average) specific star formation rate (sSFR), this 
suggests that the sSFR enhancement due to a recent major 
merger is modest. Figure [7] shows that the median sSFR of 
the bluest 30% (in the NUV - r colour) of disturbed SGs 
is around 0.15 dex (~40%) higher than the corresponding 
median value in their relaxed counterparts. A similarly mod- 
est enhancement in sSFR (<60%) due to major mergers ha s 
been reported in recent cosmological simulations (|CerJl20l"lh 
within similar stellar mass and redshift ranges as those stud- 
ied in this paper. Note that the SFR is an average value over 
the last 10 7 yr. Major mergers thus appear relatively insignif- 
icant, both in terms of driving the buildup of the SG stellar 
mass, and enhancing the star formation that is already being 
driven by other processes e.g. cold-mode accretion and minor 
mergers. While our approach is empirical in nature, our con- 
clusions regarding the overall insignificance of major mergers 
appears consistent with bo th recent theoretical work (e.g. 
ICerfalll ; iDekel et al.| [2009l and an emerging observational 
literature that indicate s a high fraction of systems with 
disk- l ike properties (e.g. [Shapiro et al]|2008l : IStockton et al.l 
|200S| ;I Forster Schreiber et al.ll2006l ) and a remarkably mod- 
est incidence of ma jor mergers (see e.g. iGenzel et al.|[200^ ; 
iTacconi et al.l l2010l . and references therein) amongst star- 
forming galaxies at z ~ 2. 

We conclude by summarising the derived star-formation 
histories of individual SGs. In the top row of Figure [8] we 
plot the starburst ages (T) and timescales (r) of individual 
SGs against their stellar masses. In the bottom row in this 
figure we show the distributions of metallicities and internal 
dust extinctions for our SGs. The yellow shaded region in 
the top-left panel shows a progressive 'one-sigma' fit to the 
data. The red line indicates the mean value, while the yellow 
shaded area indicates the region that encloses 68% of the 
galaxies. Individual error bars for each SG are shown and 
the errors in the galaxy stellar masses are typically better 
than 0.3 dex. 

We find that the star-formation timescales (r) in SGs 
are relatively short, with a median of ~0.6 Gyr, with uncer- 
tainties on individual timescales of ~0.3-0.4 dex. Given the 
large uncertainties, timescales as large as ~1.5 Gyr cannot 
be ruled out for most of the individual SGs in our sam- 
ple. However, the median value for the sample as a whole 
is smaller than 1 Gyr, consistent with previous measure- 
ments of star form ation timescales in in termediate and high- 
redshift SGs (e.g. iFerreras et al.ll2004r ) and t he high values 
of alp ha enhancement found in local SGs (e.g lThomas et all 
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In the mass range 10 Jb M© < M* < 10 iu ' 5 Mq, the star- 
burst timescales (r) and ages (T) do not correlate strongly 
with galaxy stellar mass. However, the timescales scatter 
towards lower values (< 0.3 Gyr) for more massive galax- 
10 ' 5 Mq) and a trend towards increasing age 
we move to systems with higher stellar 
masses. For example, galaxies at the upper end of our mass 
range (M« > 10 11,5 M Q ) are typically ~2 Gyrs older than 
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Figure 9. TOP: Histogram of formation redshifts, calculated 
from the starburst ages (T) and the observed redshifts of the 
SG population (see text in Section 4 for details). Median values 
are shown using the dotted vertical lines. Look-back times cor- 
responding to the formation redshifts are indicated on the top 
x-axis. BOTTOM: Formation redshift vs stellar mass in individ- 
ual SGs. The error bar indicates the average uncertainty in the 
formation redshifts across the SG population. The typical errors 
in the stellar masses are better than 0.3 dex. 



their counterparts with masses ~ 1O 1O - 5 M . These trends 
are qualitatively consistent with both observational (e.g. 
iThomas et~ai1l2005l; J uncau et al.ll2005l) and theoretical ev- 
idence (e.g. iNeistein et al.ll2006l : ICattaneo et al.ll2008r i for 
'downsizing' in massive galaxies. It is worth noting, how- 
ever, that a smooth downsizing trend with galaxy mass is 
not observed, and the large scatter in the starburst ages (T) 
indicate that SGs are not a particularly coeval population. 
While such coevality has been suggested in the past by the 
optical colours of local SGs, this is likely due to the fact that 
8-10 Gyr of evolution washes out the details of the stellar 
mass assembly. One also has to entertain the possibility that 
the trends observed in local SGs, between stellar mass and 
quantities such as luminosity-weighted age, are influenced 
by late star formation. For example, an identical gas-rich 
satellite will likely create a larger mass fraction of young 
stars, and therefore a lower luminosity-weighted age and a 
larger dilution of [a/Fe], in a smaller SG. 
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The median value for the internal dust extinction of our 
SG sample is Eb-v ~ 0.2 mag, around a factor of 2 higher 
than the typical values found in local SGs (Kaviraj et al. 
2007a). Super-solar metallicities (Z > Zq) are favoured, but 
with significant uncertainties that are typical of this type of 
photometric parameter estimation. We conclude this section 
by presenting the 'formation redshifts' of the SG population 
in Figure [9] The formation redshift is calculated using the 
derived starburst age and observed redshift of the SG in 
question , i.e. it is the redshift around which the star forma- 
tion in the galaxy is likely to peak in a WMAP-7 cosmology. 
The top panel in Figure 8 indicates that the peak of the star 
formation that builds the SGs in our sample is likely in the 
redshift range 2 < z < 5, with a median value of z ~ 3. The 
bottom panel in this figure shows the formation redshifts 
plotted against the stellar masses for our SGs. The general 
trends are similar to the plot of age (T) vs stellar mass (top- 
left panel of Figure [SJ, with the larger scatter mainly due to 
the non-linearity between look-back time and redshift. 



5 SUMMARY 

We have studied ~330 newborn spheroidal galaxies (SGs) 
in the redshift range f < z < 3, to study their high-redshift 
origin and gain insight into when and how the old stellar 
populations that dominate today's Universe formed. SGs 
have been identified by visual inspection of near-infrared 
images from the WFC3 Early-Release Science programme, 
which trace the rest-frame optical wavelengths at these red- 
shifts. Ten-filter HST photometry, covering the rest-frame 
UV-optical wavelengths, has been used to study the photo- 
metric properties of the newborn SG population and empir- 
ically estimate their star formation histories. 

The rest-frame UV-optical colours of the SG population 
indicate that virtually none of these galaxies has completed 
their stellar assembly by z ~ 5. The derived star formation 
histories indicate that the stellar assembly of our SGs likely 
peaked in the redshift range 2 < z < 5, with a median value 
of z ~ 3. Given that around hal f the present-day SG popula- 
tion w as in place by z ~ f (e.g. iBell et a l. 2005; Fabc ret al.l 
l2007f ). this implies that a significant fraction of the old stars 
that dominate the local Universe are likely to have formed 
at these epochs. 

Our results show that the star formation episodes that 
built the massive SGs are relatively short and have de- 
cay timescales less than 1.5 Gyr (with a median of ~0.6 
Gyr). Starburst ages and timescales show no correlation 
with galaxy stellar mass in the mass range 10 9 ' 5 Mq < 
M* < 10 10 ' 5 Mq. However, the timescales scatter towards 
lower values (< 0.3 Gyr) for more massive galaxies (M» > 
1O 1O ' 5 M0) and a trend towards increasing age becomes ap- 
parent as we move to higher stellar masses, with galaxies 
that have M* > 10 11,5 M Q being ~2 Gyrs older than those 
with M* < 1O 1O,5 M0. However, a smooth downsizing trend 
with galaxy mass is not observed, and the large scatter in 
the starburst ages indicate that SGs are not a particularly 
coeval population. 

Around half of the SGs in the blue cloud are relaxed, 
i.e. show no morphological disturbances of a recent merger. 
At the depth of the WFC3 images employed in this study, 
tidal debris from major mergers is likely to be visible at the 



epochs probed here. Thus, relaxed SGs in the blue cloud 
are unlikely to be driving their star-formation episodes via 
major mergers, suggesting (indirectly) that they may be ex- 
periencing minor mergers, or are the remnants of the recent 
collapse of clumpy disks, in which the star formation has 
been fed by cold streams (as envisaged in recent theoretical 
work). Furthermore, those SGs that do show tidal features, 
and are therefore likely to be recent major-merger remnants, 
exhibit only modest enhancements in their specific star for- 
mation rates of ~40%. Thus, major mergers appear rela- 
tively insignificant, both in terms of driving the buildup of 
the SG stellar mass and enhancing the star formation that 
is already being driven by other processes (e.g. cold-mode 
accretion, minor mergers.) 

This study offers empirical insights into the formation 
of newborn SGs in the early Universe, using the rest-frame 
optical data that is rapidly becoming available from new 
space and ground-based surveys. We conclude this paper 
by outlining several outstanding issues that demand further 
study. The analysis of very massive galaxies (M* > 10 Mq) 
requires a larger, statistically-significant sample of objects, 
such as the complete CANDELS survey, which offers a factor 
of ~20 increase in area compared to the WFC3 ERS pro- 
gramme (but typically in fewer filters compared to the 10- 
band HST data used here). While we have presented plau- 
sible evidence for the general insignificance of major merg- 
ers in producing SGs and driving stellar mass growth at 
high redshift, the role of such mergers in the redshift range 
1 < 2 < 3 has to be quantified further, by directly compar- 
ing the rate of emergence of SGs with the observed major- 
merger rate. Again, this will benefit from better statistics 
than is available in either the WFC3 ERS sample or the ex- 
isting CANDELS data. Bulk visual classification of galaxy 
morphologies, required for the reliable identification of merg- 
ers and SGs (see e.g. Kaviraj et al. 2007a; Kartaltepe et 
al. 2010) from large datasets such as CANDELS, can be 
achieved using novel techniques such as the Galaxy Zoo 
project. 

While rest-frame UV/optical photometry provides rea- 
sonable constraints on galaxy star formation histories, they 
can be (significantly) improved using forthcoming panchro- 
matic data. Spectroscopic line indices and/or radio contin- 
uum data can put strong priors on stellar/gas-phase metal- 
licities and the current star-formation rate respectively, en- 
abling us to significantly reduce uncertainties on the derived 
star formation history parameters (c.f. Section 3). Intrigu- 
ingly, the derived star formation timescales in our SGs ap- 
pear somewhat longer than those implied by the high [a/Fc] 
ratios observed in the low-redshift Universe. However, pho- 
tometric analyses such as the one performed here typically 
yield large error bars on the derived timescales, making it 
difficult to ascertain whether this is a real effect. A system- 
atic study of the stellar/gas- phase metallicities, [a/Fe] ra- 
tios and star formation rates of massive SGs at z > 1, using 
forthcoming spectrographs like KMOS and radio continuum 
surveys using the SKA precursors (e.g. e-MERLIN), thus be- 
comes a compelling exercise. Looking further ahead, future 
morphological studies, using the JWST and the extremely 
large telescopes, will enable us to probe massive galaxies be- 
yond z > 3, bridging the gap between SGs at the epoch of 
peak star formation and their progenitors. 

In future papers, we will systematically tackle these is- 
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sues using the forthcoming datasets mentioned above, both 
to further our understanding of the emerging SG popula- 
tion and, in particular, to bring such empirical results to 
bear on our emerging theoretical models for describing the 
high-redshift Universe. 
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